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We study the fate of the so-called ©//-loop-current order that breaks both time-reversal and parity 
symmetries in a two-dimensional hot spot model with antiferromagnetically mediated interactions, 
using Fermi surfaces relevant to the phenomenology of the cuprate superconductors. We start from 
a three-band Emery model describing the hopping of holes in the CuCH plane that includes two 
hopping parameters t pp and t p d, local on-site Coulomb interactions Ud and U p and nearest-neighbor 
V p d couplings between the fermions in the copper [Cu(3d a ,2_ M 2)] and oxygen [0(2 p x ) and 0(2 p y )\ 
orbitals. By focusing on the lowest-energy band, we proceed to decouple the local interaction Ud 
of the Cu orbital in the spin channel using a Hubbard-Stratonovich transformation to arrive at 
the interacting part of the so-called spin-fermion model. We also decouple the nearest-neighbor 
interaction V p d to introduce the order parameter of the ©//-loop-current order. In this way, we 
are able to construct a consistent mean-field theory that describes the strong competition between 
the composite order parameter made of a quadrupole-density-wave and d- wave pairing fluctuations 
proposed in Efetov et al. [Nat. Phys. 9, 442 (2013)] with the ©//-loop-current order parameter 
that is argued to be relevant for explaining important aspects of the physics of the pseudogap phase 
displayed in the underdoped cuprates. 
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I. INTRODUCTION 


The physics of the pseudogap phase of cuprate su¬ 
perconductors remains one of the most enduring open 
problems of condensed matter physics. There are re¬ 
cent pervasive hints that the pseudogap phase in most 
underdoped cuprate superconductors might involve one 
or more symmetry-breaking “hidden” orders, whose pre¬ 
cise microscopic mechanisms are still elusive to this date. 
State-of-the-art experiments such as nuclear magnetic 
resonance^, pulsed-echo ultrasound experiments^, x-ray 
scattering^ - — and scanning tunneling microscopy^ per¬ 
formed in non-Lanthanum-based materials established 
the emergence of a dome-shaped short-range incommen¬ 
surate d-wav o 9 ' 10 charge-density-wave (CDW) at low hole 
doping with a modulation described by the wavevectors 
Qx = (QojO) and Q y = (0,Qo) oriented along the prin¬ 
cipal axes of the CuCb unit cell (with Qq ~ 0.255 in re¬ 
ciprocal lattice unit s 11 ' 12 ). Quite surprisingly, the peak 
of this short-range charge order dome occurs approxi¬ 
mately at the universal hole doping x — 0.12 for sev¬ 
eral compounds—, despite their differences in material- 
specific properties. This could suggest that simple, 
low-energy effective models may potentially capture the 
essence of the physics of these materials^ - —. We will 
follow this point of view in the present work. Moreover, 
by applying pressure on these systems, the charge order 
can be completely suppressed, while the pseudogap phase 
remains unaffected— This clearly indicates that such a 
CDW order emerges on top of an already-formed pseudo¬ 
gap phase, instead of being its driving force. On the other 
hand, at very high magnetic fields, the d -wave supercon¬ 
ducting phase displayed by these materials is destroyed 


and the short-range CDW turns into a long-range order. 
In this context, it plays a central role in reconstructing 
the Fermi surface of these compounds into pockets, as is 
evidenced in quantum oscillation experiment s 18 : 19 . 

In addition to these salient features taking place within 
the pseudogap phase, another form of “hidden” order 
representing potentially one of the driving forces of quan¬ 
tum criticality in these systems (that may even co¬ 
exist with CDW order and d -wave superconductivity 
at lower temperature scales T < T*) is suggested by 
a different set of equally groundbreaking experiments: 
spin-polarized neutron scatterin g 20 : 21 and Kerr-rotation 
experiment s 22 ' — indicate spontaneous breaking of both 
time-reversal and parity symmetries in this phase at tem¬ 
peratures that are reasonably close to T* over a wide 
doping range. This phase transition has thus been re¬ 
ferred to as the Kerr transition in the literature. This 
transition was given a theoretical framework in the pro¬ 
posal by Varma— (see also an interesting, alternative 
proposal put forward in Ref—) that orbital loop cur¬ 
rent order - we shall specialize in the present work to 
the so-called ©//-phase - may account for the observed 
properties in these materials, since it naturally preserves 
the translational symmetry of the lattice and, addition¬ 
ally, it leads to the breaking of the correct discrete sym¬ 
metries consistent with spin polarized neutron scatter¬ 
ing experiments— This theoretical description requires 
starting from at least a three-band model, in which one 
includes besides the usually considered copper d x 2 _ y 2 - 
orbital, also the oxygen p x and p^-orbitals of the C 11 O 2 
unit cell. Such a minimal model turns out to be essential 
to describe intra-unit-cell loop currents involving charge 
transfer between oxygen orbitals that appear in the afore- 
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mentioned ©//-phase. This theoretical proposal is phys¬ 
ically appealing but it has one potential disagreement 
with experiments: it is hard to obtain the result that the 
underlying Fermi surface gaps out at all, since the phase 
transition does not break translational symmetry near 
the hot spots (i.e. the points in momentum space where 
the Fermi surface intersects the antiferromagnetic Bril- 
louin zone boundary). Moreover, it is important to men¬ 
tion that recently a quantum critical point (QCP) was 
revealed in the cuprates at a hole-doping x cr it — 0.18 
via an analysis of the quasiparticle mass enhancement 
using quantum oscillation experiments!!. Interestingly, 
this critical point may represent approximately the ter¬ 
mination of the Kerr transition line, the charge-order- 
dome and an as-yet-unidentified third phase competing 
with the previous two orders at a doping level reasonably 
close to optimal doping. 

On the theoretical front, the hot spot model emerges as 
an interesting, minimal low-energy effective model that 
captures qualitative aspects of the physics of the high- 
T c cuprates from a weak-to-moderate coupling perspec¬ 
tive. In this respect, an important work by Metlitski and 
Sachdev!! consisted in the elegant demonstration that, 
if the energy dispersion of this model is linearized, an 
exact emergent SU(2) pseudospin symmetry relating a 
d- wave singlet superconducting (SSC) order to a d- wave 
quadrupole-density-wave (QDW) order at wavevectors 
along the Brillouin zone diagonal (±Qo,±Qo) is veri¬ 
fied at the spin-density-wave (SDW) quantum critical 
point. This degeneracy between these two orders effec¬ 
tively produces a composite order parameter (denoted by 
QDW/SSC) with both bond order and preformed pairs 
at high temperatures as shown by Efetov et al »!! and 
the properties of this state have been explored in con¬ 
nection with the physics of the cuprates using differ¬ 
ent approaches in many works!! - —. In addition to this 
fact, another emergent SU( 2) degeneracy relating two 
additional orders - a superconducting order with a fi¬ 
nite Cooper-pair center of mass momentum (the so-called 
pair-density-wave (PDW)!! - — ) and a d-wave CDW at 
the experimentally observed wavevectors Q x and Q y 
has also been recently verified in the model in the work by 
Pepin et a/— and explored further by Wang et alA 2 . This 
additional degeneracy generates another composite order 
parameter (denoted by PDW/CDW) with similar energy 
scales that also competes with the QDW/SSC order— - —. 

In the present work, we will consider the relevant sce¬ 
nario in which yet another order parameter (the ©//-loop 
current order) competes with the QDW/SSC order in an 
effective hot spot model. The purpose of this study is to 
demonstrate the possibility that, due to this competition, 
QSW/SSC is strongly affected by the ©//-order param¬ 
eter that breaks both time-reversal and parity symme¬ 
tries, but instead preserves their product. This opens 
an interesting avenue for future research and could be 
an explanation as to why the charge-order signal in the 
cuprates is always observed along the axial vectors (i.e. 
Q x and Q y ) and never along the diagonal direction. In 


order to perform this investigation, we will construct a 
novel mean-field theory by including both ©//-loop cur¬ 
rent order and the QDW/SSC composite order parameter 
in such an effective model. As will become clear shortly, 
we will confirm in this analysis the strong competition 
between ©//-loop current order and the QDW/SSC en¬ 
tangled order, with one order parameter being clearly 
always detrimental to the other. Then we proceed to 
discuss the physical implications of this strong competi¬ 
tion for the physics of the underdoped cuprates, in light 
of the recent experiments performed in these materials. 

Technically speaking, we will introduce a three-band 
model (Emery model) describing hopping of holes in the 
Cu0 2 plane which includes two hopping parameters t pp 
and t pc i, on-site LQ and U p local interactions and nearest- 
neighbor V p d couplings between the fermions in the cop¬ 
per (d x 2 _ y 2 ) and oxygen ( p x and p y ) orbitals. By fo¬ 
cussing on the lowest-energy band, we will decouple the 
local interaction Ud of the Cu orbital in the spin channel 
using a conventional Hubbard-Stratonovich transforma¬ 
tion to arrive at the interacting part of the so-called spin- 
fermion model. Then, we will follow closely the method¬ 
ology explained in full detail in the paper by Efetov et 
alAA to define the composite order parameter associated 
with the QDW/SSC fluctuations. In addition to this, we 
will also decouple the nearest-neighbor interaction V p d 
of the model to introduce the order parameter associ¬ 
ated with the ©//-loop-current order. Lastly, we will 
proceed to derive analytically and then solve numerically 
the resulting mean-field equations, which describes the 
competition between these two order parameters. 

This paper is organized as follows. In Section II, we 
define the three-band model that we will be interested in 
and we show how to decouple the interactions to obtain 
the resulting mean-field equations describing the compe¬ 
tition between the two orders. Since the interactions that 
promote QDW/SSC and ©//-loop current order turn out 
to be different, this decoupling is unambiguous. In Sec¬ 
tion III, we solve numerically the self-consistent mean- 
field equations and then we discuss our main results. Fi¬ 
nally, Section IV is devoted to our conclusions. 

II. THE THREE-BAND MODEL 

We start this section by writing down both the non¬ 
interacting and interacting Hamiltonians of the so-called 
three-band (Emery) model following Refs. 44-47 in order 
to describe the underdoped cuprates as follows 

Ho = —t p d ^2 y^(d|. a Pi+C'/2,a + H.C.) 

i,a v 

— t pp ^2 ^2 (^i+£>/2,erP i + £ '7 2 > CT + H.C.) 
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FIG. 1: (Color online) Orbital structure and the interactions 
of the three-band model in the CUO 2 unit cell. 
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This model Hamiltonian describes the fermionic motion 
on the copper [Cu(3d x 2_ y 2)] and oxygen [0(2p x ) and 
0(2p y )] orbitals that are located in the Cu 02 unit cell 
(see Fig. [Q). The quantities d\ a , dj i<7 , p\ +0 / 2a and 
Pi+o/ 2 ,a ar e, respectively, the creation and annihilation 
operators of fermions situated on the site i with spin 
a of the Cu orbital and the the creation and annihila¬ 
tion operators of fermions on the site i + u/2 (u = x,y) 
with spin a of the O orbitals. Besides, nf a and n P i+ ~ , 2 a 
correspond, respectively, to the fermionic number opera¬ 
tors for particles located on the Cu and O orbitals. The 
model also takes into account pair hopping (t p d and t pp ), 
on-site (Ud and U p ) and nearest-neighbor ( V p d ) interac¬ 
tions involving the fermions on the Cu and O orbitals. 
The parameters Ed and e p are, respectively, the Cu and 
O orbital energies and p is the chemical potential which 
controls the electronic density in the system. 

Following Abanov and Chubukov^i, we first decou¬ 
ple the Ud part of the interacting Hamiltonian in the 
spin channel using a conventional Hubbard-Stratonovich 
transformation. The resulting action becomes 


S intl d ^} 


= xj dr^dl^i • 


+ \ [ drd 2 r 


— (dr$) 2 A (V$) 2 + m a $ 2 + §((<?) 2 ) 2 

L V s 2 


(3) 


where the bosonic field </>, = ($, 4%, is the spin- 

density wave (SDW) order parameter at the antiferro¬ 
magnetic wave vector Q = (7r,7r), v s is the spin-wave 


velocity, and m a is the spin-wave bosonic mass which 
vanishes at the quantum critical point (QCP) of the the¬ 
ory. The (7° (a = x , y , z) are the usual Pauli matrices. 
Notice that in Eq. ([5]) we have partially integrated out 
the high-energy fermions in order to derive an effective 
theory S^ t [d, </>] that corresponds to the so-called spin- 
fermion model describing the coupling between the itin¬ 
erant low-energy fermionic excitations and the antiferro¬ 
magnetic SDW fluctuations. Another possibility in order 
to investigate the Emery model is to start from a more lo¬ 
calized picture by mapping the model defined in Eqs. (1) 
and (2) onto an effective three-band t—J model^L We, 
however, will not follow this latter route in the present 
work. For this reason, we would like to state clearly from 
the outset that our starting point here will be a more 
itinerant picture. 

Now we turn our attention to the V p d interaction term 
in Eq. ©. This can be rewritten as 


Z Z = -Vpd Y Y > ( 4 ) 


i,j a,a 


with the held operators on the right-hand-side of the 
above equality being 

"^ijCT = 9 \{ d i,aPi+x/2,a A d i,crPi—x/2,cr) 

A ( d i,aPi+y/ 2,ct A d j jG Pi—y/2,a )\, (5) 

•' 4 i,er = 7) [{ d i,trPi+x/2,tr ~ d l,aPi—£/2,cr') 

A ( d i y aPi+y/2,a — d^aPi—y/ 2, it)]- (6) 

As first shown by Varma^i, only the order parameters as¬ 
sociated with ApJ, A\ 3 l, and lead to states with the 
presence of stationary-loop currents on the Cu 02 plane 
and, of course, to time-reversal symmetry breaking. The 
loop-current order with order parameter defined in terms 
of is conventionally called the ©/-loop current phase, 




FIG. 2: (Color online) Loop current pattern in the Cu02 unit 
cell for the 0/- and 0//-loop current phases [panels (a) and 
(b), respectively] proposed by Varma to explain the physical 
properties of the pseudogap state in high-Ti, cuprate super¬ 
conductors. The symbols (0) and (0) represent the orien¬ 
tation of the local magnetic moments generated by the loop 
currents. 
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FIG. 3: (Color online) Representation of the Brillouin zone 
with the underlying noninteracting Fermi surface (that en¬ 
closes the blue area) which characterizes the underdoped 
cuprate superconductors. The small black circles denote the 
so-called hot spots which are defined as the intersection of 
the Fermi surface with the antiferromagnetic zone bound¬ 
ary. For instance, the hot spot labeled as 1 has a wavevec- 
tor ki = (K-,K+) in momentum space with the constraint 
K- + K+ = 7T. The wavevectors of all the other hot spots 
in the Brillouin zone are obtained by simple symmetry oper¬ 
ations. 


where Rue 1 ® 11 = V pd a complex order pa¬ 

rameter. The mean-field value of the phase 0// was de¬ 
termined in Ref*2 4 as being equal to ±7 t/ 2. In what fol¬ 
lows, we will choose for simplicity the positive value of 
0//, since it has been shown in Ref.= 4 that this choice 
minimizes the energy for the present case. 


At this point, we would like to point out that we will 
consider only the lowest energy band of the noninteract¬ 
ing Hamiltonian defined in Eq. CD. For physically mo¬ 
tivated choices of the parameters in the present model, 
the low-energy band will naturally give rise to a Fermi 
surface shown in Fig. [3] The most singular contribu¬ 
tion in this effective model will arise from the points at 
the Fermi surface (the so-called hot spots) that represents 
the intersection of this surface with the antiferromagnetic 
zone boundary. Therefore, we will restrict the analysis of 
the present model to the vicinity of these important hot 
spot points in the considerations that follow. With this 
in mind and to set up our notation, we now define the 
following 16-component fermionic spinors 


while the loop-current order with order parameter given 
in terms of _4- 3 J, and A^l are known as the ©//-loop cur¬ 
rent phase (see Fig. ED. In view of the interpretation of 
some experiments on the pseudogap phase of the cuprate 
superconductors as an evidence in favor of the ©//-loop 
current phas o 20 ' 26 , we will analyze henceforth only this 
type of order. In this way, the decoupling of the inter¬ 
acting term in Eq. 0 using a Hubbard-Stratonovich 
transformation yields the following expression 


exp \V pd dr 


i g,g' 


= / P[i?//,0//]exp<j / drJ2 -M L +Rne l0II At {3) 


[ 2V pd 


+ Rue- i&II Afl 
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( 8 ) 

are also 


spinors in the spin space a. The symbols E, A, and L 
represent independent pseudospin spaces that are gen¬ 
erated by the Pauli matrices^. By making use of the 
results in Eq. 0, 0, ©, and ® and linearizing the 
excitation spectrum of three-band model around the hot 
spots, we obtain that the total action of the system yields 


R[Px,Py, d, (j)', Tip, Rjj ] — So [PxiPyi d] + S int \d, (j)] + S int \PxiPyi d, Tip, Rjj] 

r / d T + t; p fiTf a (—*V) fda; ^ T2 x*< 9 x \ / p x (X)\ 

= / (pUX), pI(X), dft(X)) f! + f 2 (-tV) dr + tip T ly - T 2y idy Py(X) dX 

J Ut *t.. n w f-t-o dT +ti d )\ d(X)J 


:( 2 )[ 


F t — F t if) F t — W if) 


+ (V/) 2 + 4- fO?) 2 




1 V 

dX + A 


2 y^V 


d t (X)Ei^(X)dd(X) 


dX 


( 9 ) 
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where = £ P + -£-U p — /x, ^ = Ed — /x, and both time 
and space coordinates have been collected in terms of 
the variable X = (r, r). The matrices Id, Id, T la; ( y ), and 
f 2 x(y) appearing in Eq. © are diagonal in the E ® A ® L 
pseudospin space and depend on all the parameters of the 
three-band model and also on the order parameter Rjj 
for the 0//-loop current phase (see the Appendix A to 
check their definition). Here we follow RefJ^ and intro¬ 
duce the 32-component fermionic spinors in the particle- 
hole space t as 


P x = 
Py = 


1 

71 

i 

71 


( p* ) 

\MT2Vx J 1 


pt 


( P*y \ 
\icr 2 p y J 



y|(-Px> -Pi™2) t , ( 11 ) 

-pfa 2 ) t - ( 12 ) 


In addition to the fermionic fields defined above, we also 
introduce the charge-conjugated vectors as 

* = * f 73, P x = Pln, Py = PyT3, (13) 

where 73 is the usual Pauli matrix defined in the r space. 
Hence by making use of these last definitions, the action 
in Eq. m can be naturally rewritten as 


S[P X: P y , n p , R n ] 


, _ / -d T +£, p T 3 fir 3 -f 2 (-7V) f 1x73 + f 2x *9 K \ / P x (X)\ 

= / (P x (X), Py(X), $(X)) frn-f 2 (-*V) -d T +Z P T 3 f ly T 3 +f 2y id y P y {X) dX 

J \ f-t - , _ + ft m « ' ' - / \ I 


+ A 


'S?(X)Y,i${X)aH(X) 


— - — u p) dX - 

V pd 8 p 


V r 1x T 3 
1 
2 


+ ^ 2 x^x 


r^ra 


2 s‘ u s 


dX 


r(9 T </>) 2 + (V</>) 2 


TOa ^ 2 


-5 T +^r 3 / 

+ f(7 ) 2 ^ 


(14) 


In order to derive the thermodynamical properties 
of the present model, we should first integrate out the 
bosonic field in the functional integral 

Z = J exp{-S[P x , Py,^, cj>;n p , Rn\}T>[P x , Py,^, $\ 

= J exp^S[P x ,Py,'&;n p , Ru]^V[P x , P y , \H]. (15) 

I 


However, before proceeding with that, we will neglect 
from now on the spin-density-wave interaction g in the 
present model, since this coupling effectively renormal¬ 
izes to zero under the RG flow in the low-energy limits. 
As a result, the partition function of the three-band 
model may be computed in closed form giving rise to 
the low-energy effective action 


S[P X , Pyi'&'iTlp, Rn] 

, _ / - d T + £, p T 3 f 1 r 3 -f 2 (-iV) fixT 3 +f 2 x*<9x\ /P x (X)\ 

= / (P x (X), Py(X), tt(X)) f!T 3 — f 2 (—*V) ~dr + ( p T 3 f i y T 3 + f 2y id y Py(X) dX 

J \fl x 73 + ft x ia B r\ y T 3 + tl v idy -d T + Ur 3 ) \ $ W/ 


- y J [^(X)E 1 d‘4-(A)] D{X - X') [T(X')E 1 d t T(X')] dXdX' + J - ^U p ^j 


dX. 


(16) 


Here the function D{X — X') that appears as a potential 
function in the fermionic quartic interaction is the bare 
bosonic propagator. Its Fourier transform is given by 
D(oj, k) = (w 2 /r ; 2 + |k | 2 + m a )~ 1 with m a standing for 
the spin-wave boson mass that vanishes at the QCP, v s 


is the spin-wave velocity and u> denotes the Matsubara 
bosonic frequency. 

Next, we decouple the fermionic quartic term of the 
action in Eq. (ED by using a composite order parameter 
M (X, X') for both the quadrupole density wave (QDW) 
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and the d-wave singlet superconducting (SSC) orders, as 
was described in full detail in Ref.—. This is achieved 
by considering the renormalization of bosonic propagator 
D(u), k) by the fermions at the hot spots which leads 
to the appearance of the effective spin-wave propagator 


D e ff(u, k) = ( 7 |w| + |k| 2 + m a ) _1 , where 7 is naturally 
the Landau damping term. As a consequence, the low- 
energy effective action that describes the present system 
may be represented as follows 


Seff[Px: Py, V E '\n p , Rn , M\ 

, _ / -<9 t + £ p t 3 f 1 T 3 — f 2 (—*V) Ti x t 3 + t 2 xidx\ /P x (X)\ 

= / (P X (X), P V (X), $(X)) fm-faHV) -a T + c p t 3 r ly r 3 + r 2 y id y p y (x) dx 

J \ft x 73+ft x ia B t\ v r 3 + t\yid y -d T + £ d T 3 ) \nx)j 

-i J ^{X)M(X,X')^(X')dXdX' + i J J~\X- X')Ty[M{X,X')T, 1 M(X',X)Y, 1 ]dXdX' 

♦/ft-M* 


(17) 


where we have written J(X — X') = 3\ 2 D e ff(X — X') 
instead of the spin-wave propagator in order to simplify 
the notation. The order parameter M(X,X') for the 
QDW/SSC composite order is given by 


M(X, X 1 ) = b(X, X')S 3 t U q ) , (18) 

with i = ■ (19) 


Here A + and A_ are, respectively, the d -wave singlet 
superconducting (SSC) and quadrupole density wave 
(QDW) components of the order parameter defined 
above. We also point out that the matrices u T be¬ 
long to the SU (2) group^ which lead to the constraint 
|A + | 2 + |A _| 2 = 1 involving both the SSC and QDW 
sectors. Although we have constructed an effective spin- 
fermion model for the C 11 O 2 unit cell by considering only 
the Cu atoms, we point out that the QDW/SSC order pa¬ 
rameter in Eq. (TTSl) does not lead to a charge modulation 
located on the Cu orbitals. In fact, it can be show n 16 ' 28 ’ 29 
that this composite order parameter generates a charge 
modulation with a checkerboard pattern residing on the 
oxygen O sites, which is described by incommensurate 
wavevectors with respect to the lattice. 

The effective action in Eq. Cl 3 now has a quadratic 
form and the free energy of the system can be obtained as 
follows: First one has to integrate out the fermionic fields 
in the functional integral for the partition function and 
then apply the formulae TrlnG -1 = lndet(G _1 ). Fol¬ 
lowing this procedure, we determine that the free energy 


in space-time coordinates evaluates to 


F[T, n p ,R n ,M} = _ j Trln [ G -i(x,A')]dXdX' 

+ ^ J J -1 (X — X')Tr[M(X, X')’EiM(X', X)T,i]dXdX' 

+ /(§-tKK <20) 


where the matrix G 1 (X, X’) is the Fourier transform of 
G~ 1 (ie n , k). This latter function is given by 


( -i£n + £pT 3 f 1T3 — f 2 (k) fi x T 3 -f 2x fe x 

f{r 3 — f J(k) -i£ n + i P T 3 ViyT 3 -t 2y ky 

ft x r 3 - t\ x k x f| y r 3 - f\ y k y -ie n + £ d r 3 - iM(e n , k) 

( 21 ) 

The self-consistency equation for b(X, X') is derived by 
minimizing the free energy F[T, n p , i?//, M] with respect 
to this order parameter. As a consequence, we obtain the 
following equation 



1 


G-^X.X') 

v / 0 'U T 

Ea -fit 0 


dG~ 1 (X,X l ) \ 
db(X,X>) J 


E1E3 

A 



+ J~ 1 (X-X')b(X,X') 



( 22 ) 


By performing the trace operation over the space E ® 
A 0 L (g) t for the second term on the left-hand-side of 
the equation above, the order parameter b(X , X') can be 
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simply expressed as 


b(X,X') 

= lj(X-X') Tr 
= ^ J ( x ~ x ')Tr 


vls dG~HX,X') 
G(X ’ X) db(X,X>) 


G(x,x')m 3 s 3 


-ut 



where II 3 is a projector for the three-band-model space 
which is defined as 


/0 0 0\ 

n 3 = 0 0 0 . (24) 

V° 0 V 


neglecting the SSC sector of the QDW/SSC order param¬ 
eter (A+ = 0) and applying these reduction formulas to 
the matrix G~ 1 (ie n , k), we obtain, after some algebraic 
manipulations, that det[G~ 1 (ie n , k)] evaluates formally 
to 

2 2 

det[G _ 1 (ie n , k)] = v[ m \ie n , k), (30) 

1=1 m= 1 

where D\ m \ie n , k) are well-behaved functions of the 
three-band-model parameters, which are computed in de¬ 
tailed form in Appendices B and C. Thus, by inserting 
the result displayed in Eq. (l30l) into Eq. (1271) , the mean- 
field equation for b(e n , k) in terms of V^ m \ie n , k) finally 
reads 


At this point, we will make use of the ansatz b(X, X') = 
b{X - X’) [and G(X, X') = G(X - X’)} and Fourier 
transform Eq. (1231) to momentum-frequency space. As a 
result, we get the expression 


k(An,k ) = j^Ej J ( £ « -4,k-k')Trj[G 1 {ie’ n , k')] 1 

c-/ A 


dG 1 (ie( l ,k')'( dk 
X db(e' n ,k>) J( 2 ^' 


(25) 


In order to express b(e n , k) in a convenient form, we need 
to evaluate the trace that appears in the above equation. 
This problem can be circumvented by using the following 
identity 


TrjlG-^Ck ')]- 1 

I 


dG 1 (i4,k') | 

db{e' n ,k') J 
5 det [G 1 (ie( l , k')] 


det[G _ 1 (ie(j, k')] db(e' n , k') 


(26) 


b(e n ,k) = 


3A 2 T 

16 


EE 

l,m=1 £' n 


DeffjSn ~4-k- k') 

^ W «,k') 


dT>[ m) {ie' n ,k') dk 1 

db{e' n , k') ( 2 t r) 2 ‘ 1 J 

We note that, since we have set A+ = 0, the above self- 
consistency equation describes only the QDW sector of 
the fluctuations in the present system. 

As a consequence of the result in Eq. (l30l) . we de¬ 
termine after Fourier transforming the right-hand-side of 
Eq. (l20l) that the free energy of the present model has 
the following analytical form 


F[T, n p , Rn, b} 


= -T EE ln [ V i ( i£ n,k)\ 


l.m—1 £ n 


dk 

W) 2 




3A 2 


Then, by substituting Eq. (l26ll into Eq. ([25]) . we finally 
arrive at the self-consistency equation 


x D e ff(e n ~ £' n , k-k') 


dk' 

W) 2 


Rh _ < 


V, 


pd 


p? 


(32) 


b(e n , k) 


3A 2 r y- f D e f f {s n — e' n , k — k') 
16 ,7 det[G _ 1 (ie( l , k')] 

£ n 


9det[G 1 (*e( l , k')[ dk 
db(e' n , k') ( 2 ^’ 


(27) 


where we have set J(s n — e' n ,k — k') = 3A 2 D e ff(e n — 
e^k-k'). 

We now turn our attention to the evaluation of 
det[G - 1 (*£n, k)]. In order to do this, we will need to 
use the set of determinant formulas 


det 


A B 
C D 


det (A <g) D) 


det(A) det (D-CA^B), 
[det(i)] m [det(l))] n , 


(28) 

(29) 


where we have set the volume of the system to unity. In 
order to self-consistently determine the mean-field order 
parameter Rji, we need also minimize the free energy 
with respect to it. In this way, the self-consistency equa¬ 
tion for Rjj in turn reads 

1 d'Di™'\i£ n ,k) dk 

v[ m \i£ n , k) dR n ( 27r ) 2 ' 

(33) 

The solutions of both Eqs. (13T1) and (l33l) will be obtained 
in the next section, following a numerical procedure de¬ 
scribed in great detail in Appendices B and C. 


Rn = 


V pd T 


EE 

1,171— 1 £n 


III. MEAN-FIELD RESULTS 


where A and D are, respectively, n- and m-square ma- In order to investigate the interplay between both 0/j- 
trices and det (A) is different from zero. In this way, by loop-current (LC) and QDW orders in the present three- 
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FIG. 4: (Color online) (a) Mean-field values of Ru and b as a function of the nearest-neighbor interaction V p d in the limit of 
zero temperature for A = 20. (b) Mean-field values of Ru and b as a function of the spin-fermion coupling A in the limit of 
zero temperature for V p d = 14. Both solutions in (a) and (b) were obtained by performing numerical integration in momentum 
space of the self-consistency equations given by Eqs. (EH) and (13311 with a mesh of 320 x 320 points in the Brillouin zone. Here 
m a = 1CT 2 , 7 = 10~ 5 and the other interactions are set to t p d = 1, t pp = 0.5, U p = 3, and Ed — £ P = 3. The fermionic density 
on the O orbital is given by n p = 0.6 and the position of the hot spots is such that 6 = 0.93. 


band model, we solve numerically the mean-field equa¬ 
tions for Rjj and b. The present numerical approach 
consists in the discretization of the Brillouin zone with a 
mesh of 320 x 320 points. We also make the assumption 
that the order parameter 6(e„,k) does not depend cru¬ 
cially on the frequency and momentum. In this way, we 
will only investigate the ground state properties of the 
present model, which therefore allows us to evaluate the 
Matsubara sums that appear in the mean-field equations 
exactly. We perform this calculation by either varying 
the spin-fermion coupling A or the nearest-neighbor in¬ 
teraction Vpd between O and Cu orbitals. In addition, we 
fix all other couplings in the theory. The corresponding 
results are shown in Figs. 0])a) and (b). For physically re¬ 
alistic parameters in the present model (here we choose, 
e.g., Tn a — 10 , 7 — 10 , t p d — 1, t p p — 0.5, Up — 3, 

and £d — s p = 3), we observe in Fig. Oja) that the order 
parameter Ru grows continuously from zero to positive 
values as the interaction V p d is increased for A fixed. It 
can be very interesting at this point to make a rough es¬ 
timate of the magnetic moment associated with the loop 
currents described by Ru obtained here at mean-field 
level. From Fig. 0}a), we can estimate numerically that 
the ratio of the critical parameters is given approximately 
by (R-u/Vpd) ~ 0.2. Hence, by following the same calcu¬ 
lation procedure that was explained in detail in Ref4£, 
we may conclude that the ©//-loop-current phase in our 
present theory yields a magnetic moment per unit-cell of 
approximately Mlg ~ 0.19 hb- Quite surprisingly, this 
result agrees qualitatively with the experimental estimate 
of M exp ~ 0.05/xb — O.I^b found by Fauque et alM using 


spin-polarized neutron scattering experiments. In Fig. 
|4ja), it can also be seen that the QDW order parameter 
6 , by contrast, vanishes as the interaction V p d becomes 
stronger. Moreover, one can note in the same figure a 
narrow region where both order parameters can be finite 
for moderate V p d, indicating that the present three-band 
model could in principle accommodate a coexisting phase 
involving both time-reversal (LC order) and translational 
symmetry breaking (QDW order), but as can be inferred 
from Fig. 0Ja) this apparently occurs for somewhat fine- 
tuned interactions. 

We can also analyze the behavior of the same order pa¬ 
rameters as a function of spin-fermion coupling A, when 
we keep instead the interaction V p d fixed. The corre¬ 
sponding results are depicted in Fig. 0Kb). As a result, 
we find that the LC order parameter Ru is finite below 
a threshold of A and then is clearly suppressed when this 
interaction becomes larger. Once more, the behavior of 
the QDW order parameter b is essentially the opposite 
one, namely, it grows from zero to finite values as the 
spin-fermion interaction becomes stronger. In an analo¬ 
gous way to the previous case, there is also a very narrow 
window where both phases may coexist for moderate A 
and Vpd- Despite this, the generic behavior which can 
inferred from both figures is that the LC order appears 
to be detrimental to the QDW order and vice-versa. In 
other words, we may conclude at this point that, for a 
large majority of initial choices for the couplings V p d and 
A within the present three-band model, there is a strong 
tendency for the above two orders not to coexist, at least 
at mean-field level. 
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In order to analyze the sensitivity of the above mean- 
field results to changes in the physical parameters of the 
three-band model, we have also investigated its proper¬ 
ties with respect to varying both the spin-wave bosonic 
mass m a and the orbital-energy transfer Ed — s p . As the 
strength of m a becomes larger (which corresponds nat¬ 
urally to shorter SDW correlation lengths), we obtain a 
clear tendency for both critical interactions (i.e. A c and 
Vp d ) to increase even further in our numerical data. This 
result would of course lie beyond the regime of appli¬ 
cability of a mean-field approach to the present model 
and other complementary methods that include quan¬ 
tum fluctuation effects should be used to describe such 
a regime. In addition to this, we have also examined 
the dependence of our results with respect to changes 
in the orbital-energy transfer of the model. As a conse¬ 
quence, we were able to establish numerically that, as the 
difference Ed — £ P is reduced towards zero, the critical in¬ 
teractions A c and V£ d also display a tendency to increase 
further within the present approach. 


IV. CONCLUSIONS 

In the present work, we have performed a consis¬ 
tent mean-field calculation for the three-band (Emery) 
model relevant to the phenomenology of the underdoped 
cuprates. We have shown that a low-energy effective de¬ 
scription of this model may indeed exhibit both the 0 /j- 
loop-current order first proposed by Varma^i and the 
so-called QDW which arises from an emergent SU( 2) 
pseudospin symmetry that exists in the spin-fermion 
mode l 15 i 16 . As a result, we have obtained that the above 
two order parameters have a tendency to be detrimental 
to each other, at least at mean-field level. 

We would like to point out that the mean-field values 
of the critical interactions to obtain these two phases are 
relatively large compared with some physical parameters 
of the three-band model. This is expected to be an arti¬ 
fact of the mean-field approach and, for this reason, other 
complementary methods (such as, e.g., renormalization 
group techniques that include quantum fluctuation ef¬ 
fects) should be used in order to establish a quantitative 
agreement between the present model and the experi¬ 
mental data. It is also important to mention that there 
are other works in the literature, which analyzed three- 
band models using weak-coupling diagrammatic pertur¬ 
bative calculation ;) 49 ! 50 . They have confirmed that the 
QDW order with a d -wave form factor investigated in 
the present work turns out to be more stable than the 
experimentally observed charge order with a modulation 
along the axial directions. Here, we have shown that 
the three-band model can also accommodate a ©//-loop- 
current phase that breaks time-reversal symmetry, which 
seemingly acts against the QDW order. This suggests 
an appealing scenario where the ©//-loop-current-order 
strongly competes with the QDW order, with one order 
having a tendency to suppress the other (and vice-versa) 


in the present model. This clearly opens an interest¬ 
ing avenue for future research and may help rule out re¬ 
cent competing (and mutually exclusive) interpretations 
of the universal phenomenon of the pseudogap phase dis¬ 
played in the underdoped cuprates, in light of the many 
highly-precise experiments performed in those materials 
in the last years. 
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Appendix A: Definition of the f, matrices 

The matrices Pi, f 2 , fix, ? 2 x, fiy, and f 2y that ap¬ 
pear throughout this work are defined by linearizing the 
functions of the three-band model around the hot spots 
depicted in Fig. [3] The structure of the resulting matri¬ 
ces can then be simplified by resorting to a representation 
based on Pauli matrices defined in distinct pseudospin 
spaces^, which are denoted by E, A, and L. Technically 
speaking, the pseudospin space E connects hot spots that 
can be mapped onto each other by the antiferromagnetic 
wavevector Q = Different pairs of hot spots con¬ 

nected by the wavevector Q are mapped onto each other 
by the pseudospin space A. Lastly, the pseudospin space 
L connects orthogonal quartet of hot spots. Following 
these definitions, the matrices of the three-band model 
can be simply written as 


fi = — 2t pp cos(5 Is ® 1 a ® 1l, (Al) 

T 2 = f pp (sin(5A 3 ® L 3 - E 3 ® A 3 )id x 

- t pp (smSA 3 + E 3 ® A 3 & L 3 )id y , (A2) 

fix = 7 ie -^ A3 ® i3 + 72 e ieA 30 L 3 E 3 ® L 3 , (A3) 

fax = -i 7 ie-^ A 3 ® i 3 E 3 0 A 3 + i 72 e ieA 3 ® i 3 A 3 0 L 3 , 

(A4) 

f ly = 7 ie ^ A3 - 72 e- ieA 3 E 3 ® L 3 , (A5) 

fa w = —i 7 ie^ A 3 E 3 ® A 3 ® L 3 + I 72 e - lf?A 3 A 3 , (A 6 ) 
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where 6 = (K+ — K_)/2 and Is, 1 a, and 1 l are, re¬ 
spectively, the identity matrices in the E, A, and L pseu¬ 
dospin spaces. The parameters ip, (9, 71 , and 72 are de¬ 
fined as 


, R 11 , ($ 

l “ V = 2 i ^ la \2 

tan 9 = cot ( ^ 

2 t p d \ 2 / 


71 

72 


2i pdCOS 2 ^0 + ^-sin 2 

2 & sin 2 ( 5 )+^ cos 2 



(A7) 

(A 8 ) 

(A9) 

(A10) 


Appendix B: Evaluation of the Matsubara sums for 
the mean-field equations 


1. Quadrupole density wave (QDW) order 
parameter 


In order to compute the Matsubara sum in Eq. m, 
we will consider that the QDW order parameter does 
not depend on both the frequency and the momentum. 

In this manner, we can rewrite this equation as 

l(T) = —— V V [ DeffjSn^) dV^jlS^k) dk 
16 i v[ m) (ien, k) db (2,r 

(Bi) 

where we have not written explicitly the full dependence 
of b(T) to not clutter up the notation. 

There is a subtlety to obtain the analytic continua¬ 
tion of the effective bosonic propagator D e ff(e ni k) since 
this function depends on |ct; | which is not well-defined for 
complex numbers. To circumvent that, we make use of 
the two integral formulas 


m may be rewritten as 
b(T) 

: {~2 bi dznF ( z ) De ff(~ iz ’ k ) 


3A 2 y, 

16 ^ 


dk 

W) 2 

dm) t 


l,m— 1 

1 dh\ m) (z, k) 


1 dh\ m \z, k) 




(z,k) 


db 


5 <">( 2 ,k) 


db 


(B5) 


where we have used the result "D, ( 2 , k) = 


h[ m \z,’k)h\ rn> (Zj'k) (see the Appendix C for details). 
<m) dh^ m \z, k) 

Here hj rn (z. k) and — 1 ^ — are both polynomials 

with respect to z and the relation between their degrees 
is the following 


d m )i 


deg 


<9/i[ m) (z,k) 


db 


< deg[/i| m) (z, k)]. 


(B 6 ) 


As expected, a similar inequality holds for /i| m ^ (z, k) and 
dh[ m \z.. k) 


db 


The main consequence of the result in Eq. 


(1B6I) is that we can have a series expansion of the form 

1 dh\ m \z, k) 


1 dh\ m \z, k) 


d m ) 


(z,k) 


db 


h\ m \z, k) 


db 


N 


= E 




(B7) 


where N = |dim(E ® A ® L ® r) is equal to twelve and 
k) represent both the roots of h\ m \z,k.) (1 < n ^ 
N/2) and h[ m \z, k) (N/ 2 + 1 ^ n ^ N ). The coefficients 
A ; (, ^(k) are calculated as 


1 1 IUJ 

dx 

M =- 

/ —— 

7 T 

J- oo X-lUJ 

i 

r°° dx 

sgn(w) =-i 

7 T J 

-oo X-IU 


(B2) 

(B3) 


As a result, the analytic continuation of D e //(£„,k) be¬ 
comes 


D e //(-iz,k) = —- r , ,, , 1 , , 2 , -• (B4) 

-* 7 zsgn[Im(z)] + |k | 2 + m a 

As may be easily concluded, D e //(— iz, k) is not analytic 
in the entire complex plane. Indeed, it possesses a branch 
cut (see Fig. [5]) which must be avoided when performing 
complex integration. As a result, we obtain that Eq. 


A ; ( ?(k) 


A ; ( ?(k) 


dh[ m \z : k) 
db 

dh[ m \z, k) 
dz 




*=€ i < ™ ) (k) 


dh[ m) {z, k) 
db 

dh[ m \z, k) 
dz 


^i" ) (k) 


Z= ^l!n ( k ) 



(B 8 ) 


N 


+ l^n^N. 


(B9) 


Then by substituting Eq. (IB71) into Eq. (IB5I) . we ob¬ 
tain that the mean-field equation for b(T) assumes the 
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FIG. 5: (a) Integration contour C 1 for the evaluation of the Matsubara sum appearing in the mean-field equation for the QDW 
order parameter. The crosses (x) represent the poles of the Fermi-Dirac distribution function n F (z) and the wavy-blue line at 
Im(,s) = 0 is the branch cut of D e ff(—iz, k). (6) In order to perform the Matsubara sum in this case, the integration contour 
Ci can be distorted such that it transforms into the contour C 2 that avoids the poles of n f{z) and also the branch cut. 


form 


N 


ht )= ££/A<:>(k) 

l,m—l n—1 

D e ff{-iz, k) 1 dk 


--J— $ dzn F {z ) 


27 xi 


Ci 




(217) 


2 ' 


(BIO) 


The complex integral between brackets is computed by 
changing the integration contour from Ci to C 2 (see Fig. 
0 , i.e., 


1 

27ri 


dzn F (z 


Ci 


D eff(~iz, k) 

z-cS'w 


27n 


dxn F (x ) 


-i'yx + |k| 2 +m ax +- ^ (k) 


i'yx + | k| : 2 + m a x - _ ( k ) J ’ 


(Bll) 


where = x±irj and 77 —> 0 + . At this point, we employ 
the Dirac identity 


Finally after inserting the result in Eq. (IB13D into 
Eq. (IB 101) . the mean-field equation for the QDW order 
parameter at finite temperature can be simply expressed 
as 


b(T) = 


3A; 

16 


N 


EE 

l,m =1 n =1 


[|k | 2 + m a ] 2 + 7 2 [?;(”) (k )] 2 

x ».[«,•:> wi-2 P r dI xn F (x) 

n J -00 




A& } (k) 


-00 [|k | 2 +m a ] 2 + 7 2 ® 2 

dk 




(B14) 


In the limit of T —» 0, the Fermi-Dirac distribution func¬ 
tion n F (x ) becomes the step function 9{—x). As a result, 
the integral in Eq. (IB14I) involving the Cauchy principal 
value evaluates to 


X ■ = T^S(x) + P f, (B12) 

X±U] \x J 

with V standing for the Cauchy principal value in order 
to obtain the following 


27T7 


- — f dzn F (z) 


Deff{-iz, k) 


■ rn a 


- -V 


( i-nF[6^(k)] 

[| k | 2 + ma ] 2 + 7 2 [ e W ( k )]2 Ll >" 

xn F (x) 1 


dx 


[|k | 2 + m a \ 2 + 7 2 * 2 x - £, ( ™ } (k) 


. (B13) 


lim V r dx _ XUf{x) _ - _ 

t^o J_ x [|k | 2 + m a ] 2 + 7 2 x 2 x _ ^ ; (™)(k) 


7 r 


m a 


27(|k|2 + m a )2+ 7 2[^)(k)]2 




(|k| 2 +m 0 ) 2 + 7 2 $->(k)] 


■In 


ilfS’W! 


■ rn a 


, (B15) 


where now m a = m a (T = 0) is the zero-temperature 
bosonic mass. Hence in this limit the mean-field equation 
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for the QDW order parameter is given by 


o \2 2 N 


l,m —1 n= 1 1 


x s g n fe ( Z ) ( k )] + 


(|k| 2 +m a ) 2 +7 2 fe(,™)(k)] 2 

2 7ff(k) 

^(|k| 2 +m a ) 2 +7 2 [d” ) (k)] 2 


x In 


7 lgWI 

kP +m. 


A<">(k) 


dk 

(2tt) 2 


(B16) 


where we have used the identity 6{—x) = |[1 — sgn(a:)] 
in order to simplify the above equation. 


2. (-)/ /-loop-current (LC) order parameter 


The mean-field equation for the loop-current order pa¬ 
rameter can be simplified following the same procedure 
outlined above. First of all, we transform the Matsubara 
sum in Eq. (1331) into a integral over the complex plane. 
This leads to 


Rn{T) = 


x 



dh\ m \z, k) 
dRii 


—® dzn F {z) 
2m J Ci 


•Am) 


dh[ m) {z, k) 


h[ m \z, k) 9Rii 


( z > k ) 

1 rik 

JW' 

(B17) 


Then we expand the terms between brackets in the above 
equation in a series of partial fractions, i.e., 


<9/q (m) (z,k) 


h[ m \z, k) dRii 


N 


= E 


-.(m) 

{ l,n 


(k) 


=1 2 — (k) 


1 dh\ m \z, k) 
h^iz, k) 


(B18) 


where the coefficients of the expansion Hj’^(k) are given 

by 


-(m) 

l,n 


(k) 




(k) 


dh[ m \z, k) 


dRu 

£ 

15 

uT" 

II 

dh[ m \z. k) 


dz 

*=«i ( ,» 3 ( k ) 

k > 


dRu 

*=^,» 3 ( k ) 

dh[ m \z, k) 


dz 

*=€, ( ” 5 ( k ) 


1 



(B19) 



(B20) 


Having in mind the result in Eq. (IB18D . we evaluate 
the complex integral in Eq. (IB 171) as 


1 

2m 


£ E S’( k > 


n=1 



dz 


n F (z) 


= E s £ , ( k )»de£ > (k )1 . 


n— 1 


(B21) 


Lastly the mean-field equation for the loop-current or¬ 
der parameter at finite temperature is obtained by in¬ 
serting Eq. (IB21I) into Eq. (IB17I) . Therefore this yields 


2 N 

Rum = v -f £ £ / 

I,m=ln=r k ’ 

(B22) 

In the limit of zero temperature, this equation becomes 


Rn(T = 0) = !f! £ £ / HS'tkjel-lfflk)]-^, 

I,m=ln=r ^ ' 

(B23) 

with 9(—x) being the Fermi-Dirac distribution function 
n F (x) in this case. 


Appendix C: Form of the functions X>J m ^(ie n ,k) 

In order to compute the determinant det[G~ 1 (ie n , k)] 
that appears in the main text of this work, we need 
to make use of the Eqs. (l28l) and (l29l) . Then as we 
are interested in the interplay between loop-current and 
quadrupole density wave orders, we also neglect the su¬ 
perconducting sector of the matrix u T . As a result, this 
determinant evaluates to 


2 2 

det[G -1 (i£„,, k)] = T>[ m) (ie ni k), (Cl) 

1=1 m=1 

where 'D i [ m \i£ n , k) are functions whose form will be de¬ 
termined in this appendix. 

Before proceeding with that, let us define the following 
coefficients 


ci(k x ) = V2 (^~t p d + , (C2) 

ci(ky) = \[2 t pd + i-^-ky'j , (C3) 

\/2 

02 (^:r) — 2 ( tpdkx (C4) 

02 2 ( tpdky (C5) 

which are written as a function of Cu-0 hopping t p d , the 
loop-current order parameter Rjj , and the the momen¬ 
tum distance k to the hot spots. The purpose of defining 
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these four coefficients is of course to write down the func¬ 
tions (ie n , k) in a compact form. 

We then construct the set of basis functions shown in 
Table H| from the hot-spot parameter 5 = ( K + — KJ )/2 


and the Ci(k x ), Ci(k y ) (i = 1,2). As a consequence, we 
can write explicitly the v[' n \i£ n , k) as 


J 


T>[ 1} (i£ n , k) = \{(-ie n + &)[(-*£« + Cp) 2 - tppi a i(k) + &i(k))] - P 1 (0) (k)(-*e n + - t pp P 1 (1) (k)} 

x {{-ie„ + Zd)[{~i£n + i P ) 2 - i pp (a 2 (k) + & 2 (k))] - P 2 (0) (k)(-ie n + £ p ) - t pp P 2 (1) (k)} 

- b 2 [{-ie n + £ p ) 2 - tp P (oi(k) + 6i(k))][(-ie n + £ p ) 2 - t pp (a 2 (k) + & 2 (k))]| 2 , (C6) 

V {2) (ie n . k) = \{{-ie n +^)[(-*£n + i P ) 2 - tp P ( a i (k) - b 1 (k))] - M 1 (0, (k)(-*£ n + ^ p ) - t pp M 1 (1) (k)} 

x {(-*£„ + &)[(-*£« + Cp) 2 - i pp (a 2 (k) - b 2 (k))] - M^° } (k)(-*£ n + £ p ) - t pp M^\ k)} 

- b 2 [( i£ n + £ p ) 2 - tp P (oi(k) - 6i(k))][(-i£ n + £ p ) 2 - t 2 pp {a 2 {k) - & 2 (k))]| 2 , (C7) 

v { 2 1 ] (i£ n , k) = \{(-i£ n +£ d )[(-i£ n +£p ) 2 -t 2 pp {a 3 (k)+b 3 (k))} - P 3 (0) (k)(-*£ n + Cp) -t pp P 3 (1) (k)} 

x {(-*£„ + &)[(-*£« + Cp) 2 - i pp (a4(k) + b 4 (k))] - Pj 0) (k)(—i£„ + i P ) - t pp P 4 (1) (k)} 

- b 2 [(-i£ n + i p ) 2 - *pp(a 3 (k) + 6 3 (k))][(-i£ n + £ p ) 2 - ^ pp («4(k) + & 4 (k))]| 2 , (C8) 

^ 2) (*£n,k) = !{(-*£„ +£d)[(-*£n + Cp) 2 - ^pp(«3(k) - b 3 (k))] - M^\k)(-i£ n + i P ) - t pp M^\ k)} 

X {(-*£„ + £d)[(-ie„ + Cp) 2 - t pp (a 4 (k) - &4(k))] - M^°\k)(-ie n + i P ) - t pp M^\ k)} 

- 6 2 [(-f£„ + £ p ) 2 - £ pp (a 3 (k) - & 3 (k))][(—?£„ + £ p ) 2 - t pp (a 4 (k) - & 4 (k))]| 2 , (C9) 


where we have used a second set of basis functions defined 
in Table |TT] as well as the new functions 

p/ 0) (k) = ai x (k x ) + aiy(k y ) + bi x (k x ) + bi y (k y ), (CIO) 
P; (1) ( k ) = [aj(k) + bi(k)][ai xy (k) +bi xy ( k)], (Cll) 
M/ 0) (k) = a te (fc x ) + <%(fc p ) - bi x (k x ) - bi y (k y ), (C12) 
M; (1) (k) = [a; (k) - 6i(k)][afa P (k) - h xy { k)] } (C13) 

which depend on the basis functions of both tables. 

In the main text of this paper, we have to perform 


the analytic continuation £ n —> —iz for T>\ m ' > (*£„, k). By 
observing the results in Eqs. (1C6I) (1C9I) . we conclude 
that each (i£ n , k) could be written as a product of a 
function times its complex conjugate. Therefore we can 
make the analytic continuation as follows 

V[ m \z,k) = h { ™\z,k)h { ™\z,k), (C14) 

where the functions on the right-hand side of the above 
equality are given by 


h ( i\z,k) = {(0 - £ d )[(z - i P f - *pp(oi(k) + h{k))} - Pi 0 ) (k)(2: - i P ) + t pp P 1 ( 1 ) (k)} 
x {(z - Zd)[(z - £ p ) 2 - t pp (a 2 (k) + 6 2 (k))] - P 2 0 ) (k)(2; - £ p ) + t pp P 2 ( 1 ) (k)} 

- b 2 [(z - £ p ) 2 - t 2 pp (a i(k) + &i(k))][(z - Cp) 2 - t 2 pp (a 2 ( k) + b 2 (k))], (CIS) 

h { i\z, k) = {{z + ^d)[(^ + i P ) 2 - *pp(oi(k) + Mk))] - Pl 0) (M)(z + i P ) - t pp P^\k)} 

x {{* + fd)[(z + ip? ~ t 2 pp {a 2 (k) + 6 2 (k))] - P 2 0) (k)(z + i P ) - tppP 2 (1) (k)} 

- b 2 [(z + i P ) 2 - t 2 pp (a i(k) + &i(k))][(* + i P ) 2 t 2 pp (a 2 ( k) + b 2 (k))], (C16) 

b[ 2 \z,k) = {(0 - i d )[(z - i P ) 2 - t 2 pp {a i(k) - 6i(k))] - M 1 (0) (k)(z - i p ) + t pp M 1 (1) (k)} 

x {(« - ^d)[(- + ip? ~ t 2 pp {a 2 ( k) - 6 2 (k))] - M^°\k)(z - i P ) + t pp M 2 (1) (k)} 

- b 2 [{z - i P ) 2 - t 2 pp (a r(k) - 6r(k))][( Z - ?p) 2 - t 2 pp (a 2 ( k) - b 2 (k))], 


(C17) 
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TABLE I: First set of basis functions used to represent the free energy of the three-band model. Here these functions are 
written in terms of the hot-spot parameter S = (AT+ — K-)/2 and the coefficients Ci(k x ) and c%[ky) ( i = 1,2). 


Basis function 

O'lx {kx') 

d2x ( k x ) 

&3x {kx) 

Cl4x ( k x ) 
blx{kx) 
b2x{kx) 

&3x {kx) 

&4x ( k x ) 

CLly(ky) 

Cl2y(ky) 

&3y(ky) 
CL4y{ky) 

bly(ky) 

b2y(ky) 

bs y {ky) 

b4y(ky) 
dixy (k) 

0*2xy (k) 

&3xy (k) 
d4xy (k) 
blxy (k) 

b2xy (k) 
b“3xy (k) 
b^xy (k) 


Definition 

MfclTFTMMF 

\ci(k x )\ 2 + \c 2 (k x )\ 2 
\d(k x )\ 2 + \c 2 (k x )\ 2 
\d(k x )\ 2 + \c 2 (k x )\ 2 

sin<5[|ci(fc x )| 2 - \c 2 (k x )\ 2 ] + 2 cos <5 Re [cl (^^(fe)] 
sin5[|ci(fc !C )j 2 - \c 2 (k x )\ 2 ] - 2 cosdRe[ci(fc :c )c 2 (fe E )j 
sin<5[jc 2 (A:^)j 2 - \ci(k x )\ 2 } + 2 cos 5 Re[ct(fc j; )c 2 (fc a; )j 
sin 5 [jc 2 (fe a; )j 2 - jci (fc^) j 2 ] - 2 cos S Reje^ (k x )c 2 (k x )] 

|ci(fc y )| 2 + \c 2 (k y )\ 2 

j Cl {ky ) | 2 + \c 2 (ky )\ 2 

|Cl(fcy )| 2 + \C 2 (ky)\ 2 
j Cl {ky ) | 2 + \c 2 (ky )\ 2 

sin<5[|C 2 (fey)| 2 - \ci(k y )\ 2 ] + 2cosSRe[ct(ky)c 2 (k y )} 
sin<5 [|c 2(fci/)| 2 ^ \ciik y )\ 2 ] — 2cos8Re[ct(ky)c 2 (k y )] 
siin5[jci(fcy)j 2 - jc 2 (fc !/ )j 2 ] — 2cos5Re[c|(fc y )c2(fc; / )] 
sim5[jci(fcj / )| 2 - jc 2 (fey)| 2 ] + 2cos5Re[cJ(fc !/ )c2(fcy)] 

2cos5Re[cJ (k x )ci(k y ) + c 2 (k x )c 2 (k y )\ + 2smSRe[c)(k x )c 2 (k y ) — c 2 (k x )a(k y )\ 
2cos5Re[c^ (k x )a(k y ) + c 2 (k x )c 2 (k y )\ - 2 sin 5 Re[ci(fc a: )c 2 (fc !/ ) — c 2 (k x )ci(k y )] 
2cosSRe[ct(k x )ci(ky) - c 2 (k x )c 2 (k y )] + 2 sin 5 Re[c]‘(fe ;c )c 2 (fc y ) + c 2 {k x )ci(k y )] 
2cos5Re[cJ (k x )ci(k y ) — c 2 (k x )c 2 (k y )] — 2 sin S Re [c[[( k x )c 2 (k y ) + c^(k x )ci(k y )] 
2 Re[c 2 (k x )ci(k y ) + ct(k x )c 2 (k y )] 

-2Re[ct{k x )ci{k y ) + c\(k x )c 2 (k y )) 

2 Re[c 2 (k x )ci(k y ) - c$(k x )c 2 (k y )\ 

—2 Re[c 2 (k x )ci(k y ) - c\(k x )c 2 (k y )\ 


TABLE II: Second set of basis functions needed to evaluate the free energy of the present three-band model. I 11 our notation, 
the indices l and l refer respectively to the functions aj(k) [and 6/(k)] and ai( k) [and 6j(k)[. 


I a; (k) fe((k) 

1 (k x + ky) 2 + sin 2 5(k y — k x + 2 cot 8) 2 2 sin 8[(k y + cot <5) 2 — ( k x — cot 5) 2 ] 

2 (k x + ky) 2 + sin 2 S(k y — k x — 2 cot S) 2 2 sin S[(k y — cot 8) 2 — (k x + cot 8 ) 2 ] 

3 (k x — ky) 2 + sin 2 8(k x + k y + 2 cot 8) 2 2 sin S[(k x + cot 8) 2 — (k y + cot 8) 2 ] 

4 (k x — ky) 2 + sin 2 S(k x + k y — 2 cot 8) 2 2 sin 8[(k x — cot 8) 2 — (k y — cot 8) 2 } 

1 — sin S(k x — k y — 2 cot 8) k x + k y 

2 sin 8(k x — k y + 2 cot 8) —(k x + k y ) 

3 sin 8(k x + k y + 2 cot 8) k x — k y 

4 — sin S(k x + k y — 2 cot 8) —(k x — k y ) 


h ( i\z,k) = {(0 + £ d )[(z + £p ) 2 - ip P (ai(k) - 61 (k))] - M[°\k)(z + £ p ) - t pp M[ 1] (k)} 
x {{z + £d)[(z + tpf ~ t 2 pp (a 2 (k) - 6 2 (k))] - M^ 0) (k)(z + £ p ) - t pp M^ 1] ( k)} 

- b 2 [{z + ^p) 2 - tl p (a r(k) - b^kMiz + £ p )* t 2 pp (a 2 ( k) - b 2 ( k))], (CIS) 

h 2 \z,k) = {(z - &)[(* - ^p) 2 - t 2 pp (a 3 ( k) + 6 3 (k))] - P 3 (0) (k)(- ~ £ P ) + VAA)} 

x {(z - ^d)[(- - £ P ) 2 - ipp(« 4 (k) + 6 4 (k))] - Pi 0 ) (k)(^ - ^p) + t pp Pi 1 ] (k)} 

- b 2 [(z - £p) 2 - 4(a 3 (k) + 63 (k))][(z - £p) 2 - ipp(a 4 (k) + b 4 (k))], (C19) 

h 2 \z,k) = {(z + £ d )[(z + Zp) 2 - tpp(a 3 (k) + 6 3 (k))] - P^ 0) (k)(z + £ p ) - f pp P 3 (1) (k)} 

x {(z + ^d)[(« + ^p ) 2 - ipp(« 4 (k) + b 4 (k))] - Pi 0) (k)(2 + £ p ) - t pp P[ 1] (k)} 

- b 2 [(z + dp) 2 - t pp (a 3 (k) + 6 3 (k))][(z + ^) 2 - i 2 p (« 4 (k) + b 4 (k))], (C20) 
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h i 2) (^. k ) = (0 - U)[(z - - t 2 P p{a 3 (k) - 6 3 (k))] - M^\k)(z - £ p ) + t pp M^\k)} 

x {{* ~ &z)[0 + ip) 2 ~ t 2 pp (a 4 (k) - 6 4 (k))] - Mf\k){z - £ p ) + t pp M^\k)} 

- b 2 [(z - £ p ) 2 - t 2 pp (a 3 (k) - b 3 (k))][(z - d P f - 4(a 4 ( k) - b 4 (k))], (C21) 

K 2 \z, k) = {{z + dd)[{z + £ p ) 2 - t pp (a 3 (k) - 6 3 (k))] - M^°\k)(z + £ p ) - t pp M^\ k)} 

x {(^ + ^d)[(2 + ^p) 2 -i 2 p («4(k) - b 4 (k))] - Mi 0) (k)(^ + ^ p ) -t pp M| 1) (k)} 

- + £ P ) 2 - 4(«3(k) - & 3 (k))][(* + dp) 2 4(o 4 (k) - b 4 (k))]. (C22) 


According to the approach developed in the Appendix 1, ...,1V). As h\ m \z, k) and hj m \z, k) are both sixth- 
13 for solving the LC and QDW mean-field equations, order polynomials in the variable z, their roots will be 

we need to determine first the roots of h\ m \z, k) and determined by means of numerical methods, 

k) which are denoted here as k) (n = 


* Electronic address: hermann_freire@ufg.br 

1 T. Wu, H. Mayaffre, S. Kramer, M. Horvatic, C. Berthier, 
W. N. Hardy, R. Liang, D. A. Bonn, and M.-H. Julien, 
Nature 477, 191 (2011). 

2 T. Wu, H. Mayaffre, S. Kramer, M. Horvatic, C. Berthier, 
P. L. Kuhns, A. P. Reyes, R. Liang, W. N. Hardy, D. A. 
Bonn, and M.-H. Julien, Nat. Commun. 4, 2113 (2013). 

3 D. LeBoeuf, S. Kramer, W. N. Hardy, R. Liang, D. A. 
Bonn, and Cyril Proust, Nature Physics 9, 79 (2013). 

4 G. Ghiringhelli, M. Le Tacon, M. Minola, S. Blanco- 
Canosa, C. Mazzoli, N. B. Brookes, G. M. De Luca, A. 
Frano, D. G. Hawthorn, F. He, T. Loew, M. Moretti Sala, 

D. C. Peets, M. Salluzzo, E. Schierle, R. Sutarto, G. A. 
Sawatzky, E. Weschke, B. Keimer, and L. Braicovich, Sci¬ 
ence 337, 821 (2012). 

5 A. J. Achkar, R. Sutarto, X. Mao, F. He, A. Frano, 
S. Blanco-Canosa, M. Le Tacon, G. Ghiringhelli, L. 
Braicovich, M. Minola, M. Moretti Sala, C. Mazzoli, R. 
Liang, D. A. Bonn, W. N. Hardy, B. Keimer, G. A. 
Sawatzky, and D. G. Hawthorn, Phys. Rev. Lett. 109, 
167001 (2012). 

6 J. Chang, E. Blackburn, A. T. Holmes, N. B. Christensen, 
J. Larsen, J. Mesot, R. Liang, D. A. Bonn, W. N. Hardy, 
A. Watenphul, M. v. Zimmermann, E. M. Forgan, and S. 
M. Hayden, Nat. Phys. 8, 871 (2012). 

7 J. E. Hoffman, E. W. Hudson, K. M. Lang, V. Madhavan, 
H. Eisaki, S. Uchida, and J. C. Davis, Science 295, 466 
( 2002 ). 

8 M. Vershinin, S. Misra, S. Ono, Y. Abe, Y. Ando, and A. 
Yazdani, Science 303, 1995 (2004). 

9 R. Comin, R. Sutarto, F. He, E. da Silva Neto, L. Chau- 
viere, A. Frano, R. Liang, W. N. Hardy, D. Bonn, Y. 
Yoshida, H. Eisaki, J. E. Hoffman, B. Keimer, G. A. 
Sawatzky, and A. Damascelli, Nature Materials 14, 796 
(2015). 

10 K. Fujita, M. H. Hamidian, S. D. Edkins, C. K. Kim, Y. 
Kohsaka, M. Azuma, M. Takano, H. Takagi, H. Eisaki, S.- 
i. Uchida, A. Allais, M. J. Lawler, E.-A. Kim, S. Sachdev, 
and J. C. S. Davis, Proc. Natl. Acad. Sci. Ill, E3026 
(2014). 

11 R. Comin, A. Frano, M. M. Yee, Y. Yoshida, H. Eisaki, 

E. Schierle, E. Weschke, R. Sutarto, F. He, A. Soumya- 


narayanan, Y. He, M. Le Tacon, I. S. Elfimov, J. E. Hoff¬ 
man, G. A. Sawatzky, B. Keimer, and A. Damascelli, Sci¬ 
ence 343, 390 (2014). 

12 E. H. da Silva Neto, P. Aynajian, A. Frano, R. Comin, E. 
Schierle, E. Weschke, A. Gyenis, J. Wen, J. Schneeloch, Z. 
Xu, S. Ono, G. Gu, M. Le Tacon, and A. Yazdani, Science 
343, 393 (2014). 

13 O. Cyr-Choiniere, D. LeBoeuf, S. Badoux, S. Dufour- 
Beausejour, D. A. Bonn, W. N. Hardy, R. Liang, N. 
Doiron-Leyraud, and Louis Taillefer, arXiv: 1503.02033 
(2015). 

14 A. Abanov and A. V. Chubukov, Phys. Rev. Lett. 84, 5608 
(2000); Ar. Abanov, A. V. Chubukov, and J. Schmalian, 
Adv. Phys. 52, 119 (2003). 

15 M. A. Metlitski and S. Sachdev, Phys. Rev. B 82, 075128 

( 2010 ). 

16 K. B. Efetov, H. Meier, and C. Pepin, Nature Physics 9, 
442 (2013). 

17 S. Sachdev and R. La Placa, Phys. Rev. Lett. Ill, 027202 
(2013). 

18 S. E. Sebastian, N. Harrison, and G. G. Lonzarich, Reports 
on Progress in Physics 75, 102501 (2012). 

19 N. Doiron-Leyraud, C. Proust, D. LeBoeuf, J. Levallois, 
J.-B. Bonnemaison, R. Liang, D. A. Bonn, W. N. Hardy, 
and L. Taillefer, Nature 447, 565 (2007). 

20 Y. Li, V. Baledent, N. Barisic, Y. Cho, B. Fauque, Y. Sidis, 
G. Yu, X. Zhao, P. Bourges, and M. Greven, Nature 455, 
372 (2008). 

21 L. Mangin-Thro, Y. Sidis, A. Wildes, and P. Bourges, 
arXiv: 1501.04919. 

22 J. Xia, E. Schemm, G. Deutscher, S. A. Kivelson, D. A. 
Bonn, W. N. Hardy, R. Liang, W. Siemons, G. Koster, M. 
M. Fejer, and A. Kapitulnik, Phys. Rev. Lett. 100, 127002 
(2008). 

23 H. Karapetyan, J. Xia, M. Hucker, G. D. Gu, J. M. Tran- 
quada, M. M. Fejer, and A. Kapitulnik, Phys. Rev. Lett. 
112, 047003 (2014). 

24 C. M. Varma, Phys. Rev. B 73, 155113 (2006). 

25 S. Bulut, A. P. Kampf, and W. A. Atkinson, 
arXiv: 1503.08896. 

26 B. Fauque, Y. Sidis, V. Hinkov, S. Pailhes, C. T. Lin, 
X. Chaud, and P. Bourges, Phys. Rev. Lett. 96, 197001 




16 


(2006). 

2 ' B. J. Ramshaw, S. E. Sebastian, R. D. McDonald, James 
Day, B. S. Tan, Z. Zhu, J. B. Betts, Ruixing Liang, D. A. 
Bonn, W. N. Hardy, and N. Harrison, Science 348, 317 
(2015). 

28 H. Meier, C. Pepin, M. Einenkel, and K. B. Efetov, Phys. 
Rev. B 89, 195115 (2014). 

29 M. Einenkel, H. Meier, C. Pepin, and K. B. Efetov, Phys. 
Rev. B 90, 054511 (2014). 

30 D. Chowdhury and S. Sachdev, Phys. Rev. B 90, 134516 
(2014). 

31 J. D. Sau and S. Sachdev, Phys. Rev. B 89, 075129 (2014). 

32 Y. Wang and A. Chubukov, Phys. Rev. B 90, 035149 
(2014). 

33 A. Allais, J. Bauer and S. Sachdev, Phys. Rev. B 90, 
155114 (2014). 

34 A. Allais, J. Bauer and S. Sachdev, Indian .Journal of 
Physics 88, 905 (2014). 

35 A. M. Tsvelik and A. V. Chubukov, Phys. Rev. B 89, 
184515 (2014). 

36 V. S. de Carvalho and H. Freire, Annals of Physics 348, 
32 (2014). 

37 V. S. de Carvalho and H. Freire, Nuclear Physics B 875, 
738 (2013). 

38 P. A. Lee, Phys. Rev. X 4, 031017 (2014). 

39 D. F. Agterberg, D. S. Melchert, and M. K. Kashyap, Phys. 


Rev. B 91, 054502 (2015). 

40 E. Fradkin, S. A. Kivelson, and J. M. Tranquada, Reviews 
of Modern Physics 87, 457 (2015). 

41 C. Pepin, V. S. de Carvalho, T. Kloss, and X. Montiel, 
Phys. Rev. B 90, 195207 (2014); T. Kloss, X. Montiel, and 
C. Pepin, Phys. Rev. B 91, 205124 (2015). 

42 Y. Wang, D. F. Agterberg, and A. Chubukov, Phys. Rev. 
B 91. 115103 (2015). 

43 H. Freire, V. S. de Carvalho, and C. Pepin, Phys. Rev. B 
92, 045132 (2015). 

44 V. J. Emery, Phys. Rev. Lett. 58, 2794 (1987). 

45 C. M. Varma, S. Schmitt-Rink, and E. Abrahams, Solid 
State Commun. 62, 681 (1987). 

46 W. Hanke, M.L. Kiesel, M. Aichhorn, S. Brehm, and E. 
Arrigoni, European Physical Journal Special Topics 188, 
15 (2010). 

47 M. H. Fischer and E.-A. Kim, Phys. Rev. B 84, 144502 

( 2011 ). 

48 R. Thomale and M. Greiter, Phys. Rev. B 77, 094511 
(2008). 

49 S. Bulut, W. A. Atkinson, and A. P. Kampf, Phys. Rev. B 
88, 155132 (2013). 

50 A. Thomson and S. Sachdev, Phys. Rev. B 91, 115142 
(2015). 


